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Abstract. We present high time-resolution optical energy-resolved photometry of the eclipsing cataclysmic variable 
V2301 Oph made using the ESA S-Cam detector, an array of photon counting super-conducting tunnel junction (STJ) devices 
with intrinsic energy resolution. Three eclipses were observed, revealing considerable variation in the eclipse shape, particularly 
during ingress. The eclipse shape is shown to be understood in terms of AM Her accretion via a bright stream, with very little 
contribution from the white dwarf photosphere and/or hotspot. About two thirds of the eclipsed light arises in the threading 
region. Variation in the extent of the threading region can account for most of the variations observed between cycles. Spectral 
fits to the data reveal a 10,000 K blackbody continuum with strong, time- varying emission lines of hydrogen and helium. This is 
the first time that stellar emission lines have been detected in the optical band using a non-dispersive photon-counting system. 
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1. Introduction 

V2301 Oph is an interacting binary system with an orbital pe- 
riod of 113 min (Silber et al. 1994). The spectropolarimetric 
study of Ferrario et al. (1995), made during a low accretion 
state, established that its accreting white dwarf has a magnetic 
field strength of B ~ 7 MG. On this basis, V2301 Oph is a 
polar (a magnetic cataclysmic variable, CV), with the weakest 
magnetic field of any member of this class. In these systems, 
matter flows from the donor star along a ballistic stream, un- 
til it reaches a 'threading region' where it becomes attached to 
the magnetic field of the white dwarf. The material then travels 
along the field lines, until it forms a column of material settling 
onto the surface of the white dwarf via a shock (e.g. Ferrario 
& Wehrse 1999). In some systems (eg HU Aqr, Bridge et al. 
2002) the details of the accretion geometry can be reliably re- 
lated to the shape of the eclipse lightcurve. The highly variable 
nature of the eclipse profiles in V2301 Oph, however, make the 
interpretation of the lightcurve more challenging. 

With its unique combination of high time-resolution 
and energy discrimination, S-Cam2 - ESA's energy-sensitive 
photon-counting camera - has proven particularly useful as 
a tool for studying the rapidly varying optical emission of 
CVs. Technical descriptions of the instrument may be found in 
Rando et al. (1998, 2000), while the data acquisition and reduc- 
tion methods are outlined in Perryman et al. (1999, 2001). In 
our previous studies of CVs (Perryman et al. 2001, Bridge et al. 
2002, 2003, Steeghs et al. 2003) our main tool for understand- 



No 


Date 


Cycle 


Obs 


Phase Range 




(2000) 




(UTC) 




1 


April 27 


45769—45770 


03:52-04:28 


0.802-1.115 


2 


April 28 


45781—45782 


02:17-03:03 


0.706-1.106 


3 


April 29 


45793—45794 


01:03-01:43 


0.800-1.154 



Send offprint requests to: M. A. C. Perryman: mperryma@rssd.esa.int 



Table 1. Observations of V2301 Oph obtained using S-Cam2. 
The ephemeris of Barwig et al.(1994) has been used to deter- 
mine the binary orbital phases. 



ing the spectral changes has been energy-selected lightcurves 
and their associated colour ratios. Here, we present S-Cam2 
data of V2301 Oph and relate the light curves to previous ob- 
servations and demonstrate that S-Cam2 can detect and quan- 
tify contributions from individual emission lines of hydrogen 
and helium. 



2. Observations 

The S-Cam2 campaigns were performed during April 2000 
at the 4.2-m William Herschel Telescope on La Palma, with 
the instrument mounted at the Nasmyth focus: the seeing was 
variable, ranging from good to poor. The observations were 
planned so as to include three complete eclipses, with represen- 
tative data either side. Due to constraints from other program 
targets, it was not possible to observe an entire cycle. The ob- 
servations are logged in Table 1. We estimate that V2301 Oph 
was at V ~15.7, indicating it was in a high accretion state. 
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Contact 
Points 


Phase 


Eclipse 1 


Eclipse 2 


Eclipse 3 


CI 


(0.967) 


51661.17442 


51662.11582 


51663.05722 


C2,C3 


(0.973) 


51661.17489 


51662.11629 


51663.05769 


C4 


(0.995) 


51661.17661 


51662.11802 


51663.05942 


C5 


(1.029) 


51661.17929 


51662.12069 


51663.06209 


C6,C7 


(1.035) 


51661.17975 


51662.12116 


51663.06255 


C8 


(1.065) 


51661.18211 


51662.12350 


51663.06491 



Table 2. The times for the contact points for the three eclipses. 
The times are in MJD. The error on each time estimate is ~3 
sec. 



4. Spectral properties 

The intrinsic energy sensitivity of the S-Cam2 detector per- 
mits high time-resolution spectrophotometry of astronomical 
sources. The resulting data can be analysed in terms of en- 
ergy (colour) bands and their associated ratios (Perryman et 
al. 2001) or even direct spectral fitting (de Bruijne et al. 2002, 
Reynolds et al. 2003). In this paper we extend this spectral- 
fitting to accommodate the presence of emission lines, a par- 
ticularly useful technique for cataclysmic variables due to the 
swiftness of the eclipse transitions. First, however, we exam- 
ine the data in terms of energy bands, since this is useful for 
characterising broad changes in emission properties. 



In Fig n we show 'white-light' lightcurves for the three 
observations, after barycentric correction, background subtrac- 
tion, phase-folding and correction for atmospheric extinction; 
see Perryman et al. (2001) for a description of these data re- 
duction steps. The data have also been corrected for varia- 
tions in pixel-to-pixel detection efficiency and energy response. 
The light curves show that the accretion stream in V2301 Oph 
is particularly prominent in the optical band, compared with 
many other polars. 



3. Light Curves 

Our data were phased with respect to the orbital ephemeris 
of Barwig et al. (1994), which was also the ephemeris used 
by Steiman-Cameron & Imamura (1999). There are significant 
variations between the 3 eclipse profiles in each of the eclipses. 
They do, however, exhibit a similar shape to those seen in other 
studies. The first is similar to the 1997 June 2 eclipse seen 
by Steiman-Cameron & Imamura (1999), with a pronounced 
two-stepped descent into the eclipse. This two-stepped profile 
is also seen in the UV data of Schmidt & Stockman (2001). 

Steiman-Cameron & Imamura (1999) identify eight con- 
tact points in the lightcurve, associated with the ingress/egress 
of the white dwarf photosphere and the ingress/egress of the 
accretion stream. By reference to Fig 6 in Steiman-Cameron & 
Imamura, we determined the phases of their contact points and 
then plotted these phase points as vertical dashed lines in our 
Fig 1. They associate Cj, C2, C5 and C6 with the eclipse of the 
photosphere, while they associate the four other contact points 
with the eclipse of the accretion stream. We show the observed 
times of these contact points in Table |2 

The contact points as defined by Steiman-Cameron & 
Imamura (1999) generally match up with the same feature in 
our light curves. For contact points C4,g the features line up very 
well. There is some discrepancy in the eclipse ingress because 
the accretion stream in V2301 Oph is relatively bright com- 
pared to many other polars. In the second eclipse, the stream 
accounts for roughly 3/4 of the optical emission. We believe 
that discrepancies between the other contact points can be ex- 
plained by a variation in the amount of material located in dif- 
ferent regions of the accretion stream trajectory. We discuss this 
in greater depth in §5. 



4.1. Colour bands 

The data were segregated into three bands corresponding to 
blue (3200 — 4700 A), visual (4700 — 5500 A) and red (5500 
— 7200 A colours). We label these bands B', V and R', em- 
phasising that they are not on any formal photometric system. 
The boundaries between the bands are selected to ensure that 
roughly equal numbers of events are assigned to each band. In 
Fig 2, we again show the three eclipses, but this time split into 
these three energy-selected lightcurves. The amplitudes of the 
V and B' band curves are adjusted slightly to correct for small 
differences in total counts (the initial segregation only giving 
approximate equipartition among events) and the V and B' 
bands are then offset vertically on the same plots. Differences 
in slope between the bands are therefore meaningful. Although 
there are no dramatic differences between the bands, in all three 
eclipses the R' band data exhibits noticeably more structure, 
with the steps and the ledges between them more sharply de- 
fined than in the other bands. In particular, the first and third 
eclipses both contain dips in the R' band at phase ~ 0.0375, 
which are much less apparent in the B' and V data. In the first 
eclipse, as well, the ledge between the two ingress steps is more 
or less horizontal, whereas it shows a pronounced slope in the 
B' and V data. This suggests that if the material responsible 
for the steps is concentrated in two distinct regions, then the 
region between them is emitting mainly at shorter wavelengths 
than that seen in the R' band data, and/or is preferentially ab- 
sorbed in the red. 

The relatively small numbers of events per time bin in each 
colour band make the study of colour ratios difficult. In order 
to remove short-term variations due to seeing, guiding errors 
and intrinsic flickering, we first smoothed the individual R' and 
B' band lightcurves using a sliding cell with a width of 15 s. 
The ratio R'/B' against binary phase, for all three eclipses, is 
shown in the narrower panels in Fig 2. It is apparent that there 
are significant differences in the variation of R'/B' between the 
three observations, particularly in the late stages of the ingress, 
just before the onset of deep eclipse. 

4.2. Spectral fitting 

We now examine the same data using the complimentary ap- 
proach of model fitting to time-resolved spectra. Each com- 
plete eclipse observation was sliced into consecutive 25 s 
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Fig. 1. Phase-folded light curves of the three V2301 Oph observations. The right-hand plots are enlargements of the eclipse 
intervals in the left-hand plots. These white-light (A ~ 3200 - 7200 A), background-subtracted light curves have been corrected 
for variations in pixel-to-pixel gain and responsivity. A time (i.e., zenith angle) dependent correction for atmospheric extinction 
has also been applied. The data have been binned into 1 s intervals. The phases and cycle numbers are with respect to the 
ephemeris of Barwig et al. (1994). The vertical dashed lines indicate the phases of the contact points identified by Steiman- 
Cameron & Imamura (1999); ci and C5 are the contact points of photospheric ingress and egress respectively. 



intervals, with a single spectrum created from each inter- 
val. The duration of the intervals was chosen as a compro- 
mise between time-resolution and S/N, and is short enough to 
probe spectral variations on the typical timescale of the ma- 
jor ingress/egress features seen in the lightcurves. The spectra 
were corrected for electronic artifacts and pixel-to-pixel gain 
variations. Corresponding background spectra were obtained 
by using 150 s data intervals from the deepest part of the eclipse 
in all three observations, ensuring that any residual contribution 
from the donor star is removed from the spectrum, in addition 
to the sky. 

Spectral fitting made use of the XSPEC model-fitting pack- 
age (Arnaud 1996), developed for use in high-energy astro- 
physics, but adaptable to the analysis of low-resolving power 
optical data such as S-Cam spectra. In the left-hand panels of 
Fig 3, we show the result of fitting a blackbody model against 
a post-eclipse slice (the final 25 s interval) from the first obser- 
vation, together with the associated residuals. The optical flux 
therefore originates from the white dwarf, the accretion region 
and the accretion stream. The reduced x 2 f° r this fit was 4.8. 
The systematic structure in the residuals clearly indicates the 
inadequacy of this single-component model. 

In the right-hand panels of Fig 3, we show the improve- 
ment to the fit by the addition of emission lines to the black- 



body model. We fixed the blackbody temperature to be 10000 
K which was the mean value derived after experimenting with 
kT hh as a free parameter (the range in kT bb was a few 1000K). 
The lines included in our fit were Ha and H/3, both of which are 
seen in emission in the data of Barwig et al. (1994), Silber et al. 
(1994), Ferrario et al. (1995) and Simic et al. (1998), and Hei 
5875, observed strongly in emission in many AM Her spectra, 
including V2301 Oph (Barwig et al. 1994). Hen 4686 is also 
likely to be a strong feature of the spectrum. This line, however, 
lies sufficiently close to H/3 that, with our limited resolution, 
the two lines are blended. While our fit statistics show that it 
is necessary to have an emission component in this part of the 
spectrum, we are more than likely modelling the contribution 
of two lines with our single component labeled 'H/T. 

Using the published optical data of Simic et al., we esti- 
mated that the width (cr) of Ha is ~ 170 A. The line widths 
in our model were therefore fixed at this value, and the cen- 
troid wavelengths were fixed at their laboratory values. We also 
added a 'Broad Emission Component' centered on 3540 A. 
This is to take into account the many features expected to exist 
to the blue side of He n 4686, including the Balmer limit, which 
it would be impractical to fit as individual components given 
our limited resolution. The width of this component was al- 
lowed to vary, but for the other lines the only adjustable param- 
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Fig. 2. Energy-band selected lightcurves for the three eclipses 
in Fig 1, with the ratio R'/B' against phase plotted below each 
observation. The vertical dotted lines in the R'/B' plots indicate 
the approximate start and end phases of the deep eclipse. Not 
only is the eclipse ingress and egress shape in general more 
sharply defined in the R' band data than the B' band, but there 
are significant differences between the three ratio plots. 



eter was the normalisation. With all model parameters present, 
the mean reduced^ 2 over the first eclipse observation was 1 .67, 
compared to 6.0 for the blackbody alone. 

In Figure @] we show the variations in line strengths with 
binary phase for the first eclipse observations. The mean^ 2 for 
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Fig. 3. Left: blackbody fit, with associated residuals, for the fi- 
nal 25 s slice in the first observation. Right: the same data, but 
now modelled with a blackbody and emission lines. The resid- 
uals are plotted on the same scale as in the left-hand box. 

the second and third observations is worse than for the first: 
increasing to 1.9 and 2.4 respectively. This is probably due to 
long-term drifts in the gain, after warming and cooling of the 
cryogenic system. 

The top panel of Figure|4]shows the variation in the normal- 
isation of the blackbody, together with a smoothed version of 
the same trend. It is seen that this curve closely resembles the 
shape of the white light lightcurve: this is expected since most 
of the photons detected will come from the continuum rather 
than the lines. 

The trends in the three narrow lines in the model are plot- 
ted below the blackbody normalisation on the same time axis. 
In order to make the line strengths meaningful, we normalised 
them to the blackbody continuum strength: the plotted inten- 
sities represent the height of the line above the continuum, as 
a fraction of the continuum height at 4130 A. If the line in- 
tensities are proportional to the underlying blackbody normal- 
isation therefore, the resultant trend will be approximately flat. 
This is seen to be the case for Hei 5875 and H/3. Ha, how- 
ever, becomes progressively fainter relative to the above lines, 
almost vanishing around 0=0.95. It is likely that H/3 also be- 
comes weaker, but that this is masked by the contribution from 
Hen 4686 (see above). One physical reason for the Balmer 
lines to become weaker as the phase nears eclipse is that they 
are more likely to originate from the heated face of the sec- 
ondary star (and hence not visible at eclipse), while the Hen 
4686 line originates from the accretion stream, which is visible 
for longer. This also explains why Ho- is also somewhat slow 
to reappear after the sharp rise in continuum intensity during 
the egress, displaying a ~ 0.02 lag in phase. The final model 
parameter, the Broad Emission Component at 3540 A varies 
erratically, but is generally stronger in the post-eclipse phases 
than prior to the eclipse. 

5. Discussion 

In §3 we showed that the optical emission from V2301 Oph 
has a larger contribution from the accretion stream compared 
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Fig. 4. Continuum and line strengths in observation 1 . From top 
to bottom: phase-resolved, normalised line intensities for the 
following model components: blackbody (for T e g = 10000 K), 
Ha, He i 5875, and H/3. The line strengths are normalised to 
the instantaneous blackbody intensity, and refer to the height 
of the line as a fraction of the continuum height at 4130 A. The 
final model line, the Broad Emission Component, is not plotted, 
since its identification is uncertain and it varied erratically with 
time. 



to many other polars. In order to explore this further, we trace 
the Roche potential out of the binary system along the line of 
sight from any point within the vicinity of the white dwarf. We 
then adjust the system parameters, (mass ratio (q) and incli- 
nation (i), and accretion spot location), under the constraints 
of grazing contact of this line of sight with the Roche lobe at 
particular phases. We initially adopted (q,i) values of Steiman- 
Cameron & Imamura (1999), namely 0.25 and 81.75° respec- 
tively. However, with the usual assumption that the stream fol- 
lows a ballistic and then a magnetic trajectory, we found that 
the stream would be fully eclipsed by 0=0.99, whereas our 
observations show that the stream remains visible until phase 
~0.995. In order to resolve this, the donor star must become 
smaller, with a corresponding increase in the inclination. There 
are no individual constraints on (q,i), only the pairs set by the 
eclipse length, as illustrated by the curve in Fig 7. The choice 
of q=0.15, i=84.0° is sufficient to ensure stream visibility at 
phase 0.99, so in the absence of further constraints we base our 
subsequent discussion on these values. 

In Figure |6] we show the orientation of the binary at dif- 
ferent eclipse phases, for two different configurations of the 
threading region. The parameters used to generate these plots, 
apart from (q,i), are a spot azimuth of 65° and the magnetic 
field aligned at 20° to the white dwarf spin axis, both of which 
are in line with the values in Steiman-Cameron & Imamura 
(1999). The threading of material onto the magnetic field lines 
takes place between 0.2 — 0.12 a (where a is the binary sepa- 
ration) for the set of plots on the left, and between 0.4 — 0. 12 
a for the set of plots on the right. 




Inclination (degrees) 

Fig. 5. The relationship between mass ratio (q) and inclination 
for V2301 Oph, as defined by the eclipse geometry. The phase 
of the stream visibility favours q=0.15, i=84.0°. 

The scenario on the left can then be related to the eclipse 
lightcurve seen in the first observation (Fig 1). By 0=0.967, the 
white dwarf is fully eclipsed. This causes only a small dip in 
the lightcurve, however, and is only really apparent in the first 
observation. By 0=0.973, the bright part of the magnetic stream 
closest to the white dwarf has also been eclipsed. The intensity 
then declines only gradually until 0=0.985, when the threading 
region begins to be occulted. Brightness declines rapidly, and 
by 0=0.995 the stream is entirely eclipsed. 

On egress, the first component to emerge is the white 
dwarf, accounting for the sudden, steep rise in intensity. This 
is followed immediately by the hot stream, characterised by a 
slightly shallower slope. There is then another plateau as the 
empty mid-arch of the magnetic field is revealed. At 0=1.035 
the threading region again begins to emerge, and by 0=1.054 
the threading region is completely visible. The implication of 
this observation is that around one third of the stream bright- 
ness is in the region around the white dwarf, while the remain- 
ing two thirds arises in the threading region. 

Turning to the scenario on the right, an increase in the ex- 
tent of the threading region can explain the eclipse shape of the 
second observation. At 0=0.973, the threading region is already 
starting to be eclipsed, and by 0=0.985 is substantially hidden. 
However, the mid-arch of the magnetic field contains more ma- 
terial and is visible for longer. This accounts for the longer, 
brighter decline to eclipse seen in the second observation. As 
before the stream is fully eclipsed by 0=0.995. At 0=1.054, 
during the egress, the threading region is still not completely re- 
vealed, and hence the lightcurve continues to increase in bright- 
ness beyond this phase, in contrast to the first observation. 

We therefore have a self-consistent explanation for the 
white-light variations seen in the first two observations, within 
the context of the standard ballistic-plus-magnetic stream. 
Almost all of the observed light arises in the stream, and 
two thirds of that is contributed by the threading region. 
Observation three closely resembles observation one, except 
that the threading region appears somewhat less bright. 

The colour information, as presented in Section 4.1, also 
supports the above interpretation. Fig 2 shows that until the 
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Fig. 6. The eclipse geometry of V2301 Oph, forq=0.15, i=84°. 
The spot azimuth is 60°, and the magnetic field is aligned at 
20° to the white dwarf spin axis. The plots on the left show six 
phases of the eclipse for the case where the threading occurs 
later in the stream, while the plots on the right show the same 
six phases for the early threading scenario. 



stream becomes rather faint at 0=0.99, the second observation 
is redder than the first. This shows that what we are seeing in 
the second observation is cooler material in the furthest reaches 
of the threading region, whereas in the first observation, we are 
seeing material plunging deep into the field before threading. 
Observation three (in which there was a slight reddening com- 
pared to the first observation) is a little harder to understand, 
but the threading region is fainter here anyway, so is perhaps 
smaller and cooler than in the first observation. 



6. Conclusions 

We have presented energy-resolved STJ photometry of the 
eclipsing polar V2301 Oph, reporting the first detection and 
analysis of stellar emission lines based on non-spectrally dis- 
persed optical data. Three eclipses are observed in total, with 
significant variations in shape observed between them. These 
variations are shown to be consistent with the usual AM Her ac- 
cretion scenario, whereby mass flows from the donor star to the 
white dwarf along a ballistic stream, before becoming threaded 
onto the magnetic field lines of the white dwarf. By varying the 
extent of the threading region along the ballistic stream, the ma- 
jor features of the observations can be accounted for, provided 
the mass ratio is decreased to 0.15, with a corresponding in- 
crease in the inclination to i=84.0°. The colour variations seen 
in the energy-resolved data are also consistent with this inter- 
pretation, since they suggest that the material is at its hottest 
when it penetrates most deeply into the magnetic field before 
becoming threaded. 

These observations show that current STJ detectors can re- 
solve optical emission lines. However, the resolution of such 
detectors is much poorer than can be obtained with dispersive 
spectrometers. Until now, nearly all STJ work has used the ma- 
terials niobium and tantalum. The obvious next step is to mi- 
grate the technology to even lower critical-temperature super- 
conductors such as hafnium (T c = 0. 128 K, compared to 4.48 K 
for tantalum). With a band-gap well below a meV (0.021 meV), 
hafnium has the potential to provide better energy resolution at 
all energies and is therefore of keen interest, for both visible 
and X-ray applications. The predicted tunnel-limited resolving 
power of hafnium is ~64 at 1 eV, so that the wavelength reso- 
lution is ~2 A at Ly a, -50 A at 5000 A, and 150 A at 20000 
A. 
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